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We report the new measurement of D°—D mixing in decays to K + K~ 
and 7r + 7r~ final states that is based on the total Belle data sample of 
976 fb -1 . The preliminary results are yep = (1.11 ± 0.22 ± 0.11)% and 
A T = (-0.03 ± 0.20 ± 0.08)%. 
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1 Introduction 



Mixing of neutral mesons occurs when the flavor eigenstates differ from the physical 
mass eigenstates of the meson-antimeson system. In case of D° mesons mass eigen- 
states are expressed as \Dl 2 ) — p\D°) ± q\D°), with p 2 + q 2 = 1. If p = q the two 
mass eigenstates are CP-even and CP-odd; otherwise CP is violated. Charm mixing 
is characterized by two parameters: x = Am/T and y = Ar/2r, where Am and Ar 
are the mass and decay width differences of mass eigenstates, respectively, and T is 
the average D° decay width. 

Mixing in D° decays to CP eigenstates, such as D° — > K + K~, manifests in a 
lifetime that differs from the lifetime of decays to flavor eigenstates, such as D° — > 
K~tt + . The quantity 



Vcp 



t(K+K- 



is equal to the mixing parameter y if CP is conserved. If CP is violated the lifetimes of 
D° and D° 
defined as 



D° and D° decaying to the same CP eigenstate also differ and the lifetime asymmetry. 



. t(D° -> K'K+) - t(D° K+K-) 

= — =n (2) 

t(D K-K+) + r(D° K+K-) 

becomes non-zero. The quantities yep and Ar are connected to the mixing parameters 
x and y by [1] yep = y cos — ^A M xsm.cf) and A r = \A M y cos — xsin0, where 
4> = arg(g/p) and A M = 1 - \q/p\ 2 . 

First evidence for D° — D mixing was obtained in 2007 by Belle [2J and by 
BaBar [3] in two different decay modes. We report here an update of our first- 
evidence measurement in D° — > K + K~ ,tt + it~ decays using almost twice larger data 
set. 



2 Event selection 

The decays D° — > K + K~ ,tt + tt~ and D° — > K~n + are reconstructed in the decay 
chain D* + — > D°tt + in order to suppress background and to tag the D° flavor at 
production. To reject D* + candidates coming from B decays, we require that the 
D* + momentum measured in the center-of-mass system (CMS) of the e + e~ collisions 
be greater that 2.5 GeV/c; for a fraction of data taken at Y(5>S) we increase this 
threshold to 3.1 GeV/c. 

Besides using our standard kaon and pion selection criteria, based on dE/dx, 
threshold aerogel Cherenkov counters and time-of-flight, and vertex fits, we select D° 
candidates using two kinematic variables: the invariant mass of the D° decay products 
M and the energy released in the D* + decay q = (Md» — M — m n )c 2 . The proper 
decay time of the D° candidate is calculated from the projection of the vector joining 
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the two vertices, L, onto the D° momentum vector: t = m D oL ■ p/p 2 , where m D o is 
the nominal D° mass [I]. The decay time uncertainty at is evaluated event-by-event 
from the error matrices of the production and decay vertices. 

The sample of events for the lifetime measurements is selected using AM, Aq 
and at- These selection criteria are optimized on Monte Carlo (MC) simulation by 
minimizing the statistical error on yep- Background is estimated from sidebands in 
M; the sideband position has also been optimized in order to minimize the systematic 
uncertainty. The yields of selected events are 242 x 10 3 (K + K~), 114 x 10 3 (7r + 7r~) 
and 2.61 x 10 6 (i^vr+), with signal purities of 98.0%, 92.9% and 99.7%, respectively. 

3 Lifetime fit 

By studying the proper decay time distribution of D° — > K~ir + decays we observe a 
significant dependence of its mean value on cos 9*, where 9* is the D° CMS polar angle. 
Using MC simulation we find that the generated proper decay time distribution of 
selected events agrees well in each cos 9* bin with an exponential distribution and that 
the lifetime is consistent with the generated value. The observed dependence is thus 
due to the resolution function offset that depends on cos 9*. Therefore, to reduce 
systematic uncertainties arising from the resolution function parameterization the 
measurement is performed in bins of cos#*; an additional requirement | cos^*! < 0.9 
is imposed to suppress the events with the largest offsets (about 1% of events). 

The fitting procedure in each cos 9* bin is similar to the one used in our previous 
measurement [2], where we performed a binned simultaneous fit to KK, Ktc and tttt 
samples. The resolution function is constructed similarly, from normalized distribu- 
tion of a t using a double or triple Gaussian PDF for each a t bin. The widths <r£ u11 
and fractions w k of these Gaussian distributions are obtained from fits to the MC 
distribution of pulls (t — t gcn )/a t . The procedure is repeated for each cos 9* bin. The 
parameterization reads: 

"bin n 3 

R(t) = 53 h w kG{t] fa, a ik ), (3) 

i=l k=l 

where G{t;fa,aik) represents a Gaussian distribution of mean fa and width aik, and 
fi is the fraction of events in the i-th bin of the at distribution. The mean and width 
are parameterized as: 

n 

&ik = s k al ull ai fa = t + a(a % - fj a j) ( 4 ) 

where s k , k = 1, ...,n g are the width scaling factors for each of the n g Gaussian's, to 
is the resolution function offset, and a is a parameter to model a possible asymmetry 
of the resolution function. The parameters s k , to and a are free parameters in the fit. 



2 



The proper decay time distribution is parameterized with 



N r 

/W = 7 /e- 



t,/T R(t-t')dt' + B(t), 



(5) 



with the following free parameters: N, r, Sk, to an d a- A sideband subtracted at 
distribution is used to construct R(t). The term B(t) describes background and is 
fixed with a fit to the sideband distribution. 

Background is parameterized as a sum of two lifetime components, a component 
with zero lifetime and a component with an effective lifetime r&: 



The background resolution function Rb(t) is assumed to be symmetric (a = 0) and is 
composed of three Gaussian's with S3 = S2- The fraction / of zero-lifetime component 
is found to be cos 9* dependent; its value is fixed in each bin using MC simulation. The 
parameters to, s±, S2 and r& are determined from a fit to sideband distribution summed 
over cos 9* bins. However, the shape in individual bins remains bin-dependent since 
a t distribution, / and iV& depend on cos#*. 

To extract yep and A-p the decay modes are fitted simultaneously in each cos 9* 
bin using a binned maximum likelihood fit. The parameters shared between decay 
modes are yep and A? (KK and 7T7t), to and a (all decay modes) and parameters Si, 
S2 and S3 up to an overall scaling factor. The fit has been tested with the generic 
MC simulation equivalent to six times the data statistics. The fitted yep and Ar 
are found to be consistent with the input zero value and the fitted Kir lifetime r is 
found to be consistent with the generated value. A linearity test performed with MC 
simulated events re-weighted to reflect different yep shows no bias. 

4 Results 

The experimental data were taken with two different silicon vertex detector (SVD) 
configurations: for the first 153 fb _1 a 3-layer SVD was used, while for the rest of the 
data a 4-layer SVD was used. We treat both running periods separately, since the 
resolution function differs. 

The proper decay time distributions are fitted simultaneously as discussed in the 
previous section. The results of the fits are shown in Fig. [TJ Fit confidence levels 
(CL) are above 5% (except one with CL=3.3%) and are distributed uniformly El- The 
plots in Fig. [1] are obtained by summing the histograms and functions in all fitted 

*We use Pearson's definition of x 2 an d take only the bins with the fitted function greater than 
one. 



B{t) = N b f [f5(t') + (1 - /)I e -*7n ]i4( £ _ f)dt 



(6) 



n 
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cos#* bins. The residuals show no significant structure. The normalized Pearson's 
X 2 are 1.01 (SVD1) and 1.16 (SVD2). 

Fig. |2]shows the results of the fits in bins of cos 6* for yep, A-p and the D° — > K~tt + 
lifetime r. The average is obtained by a least square fit to a constant. We find 
y C p = (l.H ± 0.22)%, A r = (-0.03 ± 0.20)% and r = (408.46 ± 0.54) fs, where the 
errors are statistical only. The results for individual running periods are consistent 
with each other. The measured D° lifetime is also consistent with the current world 
average jl]. 
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Figure 1: Simultaneous fit of SVD1 (top) and SVD2 (bottom) proper decay time 
distributions. The plots show a sum of distributions and fitted functions (in red) 
over cos 9* bins; the residuals are plotted beneath. The signal-region distributions 
are shown as error bars and the sideband-region distributions as hatched histograms. 
The "(+)" and the "(-)" denote the charge of the tagging slow pion. 



5 Systematics 

The systematic uncertainties are summarized in Table [TJ The largest contribution is 
found to arise from possible SVD misalignments. The impact of misalignments has 
been extensively studied using a special signal MC simulation with different local and 
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Figure 2: Results of simultaneous fits in bins of cos6>* (points with error bars) for 
SVD1 (top) and SVD2 (bottom). The horizontal line in each plot is a least square fit 
of constant to data points. 

global SVD misalignments. We find that the SVD misalignment affects the resolution 
function considerably, especially its offset t , and it can explain the differences seen 
between MC and data. However, the effect on the resolution function is very similar 
for KK, Kn and nn decay modes, and thus mainly cancels when relative lifetime 
measurements such as yep and Ar are performed. 

Other sources of systematics arise from possible impact of selection criteria to the 
acceptance, from the position of mass window AM, from background parameteriza- 
tion and from background statistical fluctuations, from resolution function parame- 
terization and from binning in t and cos#*. The total systematic uncertainties are 
estimated to 0.11% for y CP and 0.08% for A r . 



6 Summary 

With the full Belle data set of 976 fb _1 , we obtain the preliminary results 

ycp = (+1.11 ±0.22 ±0.11)%, 
A r = (-0.03 ±0.20 ±0.08)%, 



(7) 
(8) 
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Table 1: Systematic uncertainties 



source 


Ay C P (%) AA r (%) 


acceptance 

SVD misalignments 

mass window position 

background 

resolution function 

binning 


0.050 0.044 
0.060 0.041 
0.007 0.009 
0.059 0.050 
0.030 0.002 
0.021 0.010 


sum in quadrature 


0.11 0.08 



where the first error is statistical and the second is systematical. The significance of 
Ucp 7^ is 4.5a when both errors are combined in quadrature, and 5.1a if only the 
statistical error is considered. The result for A-p is consistent with no indirect CP 
violation. Both results are in good agreement with our previous measurement [2] as 
well as with the BaBar measurements in the same decay modes [U E] . 
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